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Abstract 

This paper presents novel analytic expressions for the Rice le—function, Ie{k,x), and the incom¬ 
plete Lipschitz-Hankel Integrals (ILHIs) of the modified Bessel function of the first kind, Iem,,nici, z). 
Firstly, an exact infinite series and an accurate polynomial approximation are derived for the Ie{k,x) 
function which are valid for all values of k. Secondly, an exact closed-form expression is derived for 
the Iem,n{a, z) integrals for the case that n is an odd multiple of 1/2 and subsequently an infinite 
series and a tight polynomial approximation which are valid for all values of m and n. Analytic upper 
bounds are also derived for the corresponding truncation errors of the derived series’. Importantly, these 
bounds are expressed in closed-form and are particularly tight while they straightforwardly indicate that 
a remarkable accuracy is obtained by truncating each series after a small number of terms. Furthermore, 
the offered expressions have a convenient algebraic representation which renders them easy to handle 
both analytically and numerically. As a result, they can be considered as useful mathematical tools that 
can be efficiently utilized in applications related to the analytical performance evaluation of classical 
and modern digital communication systems over fading environments, among others. 
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Index Terms 

Rice /e—function, Incomplete Lipschitz-Hankel Integrals, Marcum Q—function, approximations, Bessel 
functions, fading channels, performance evaluation, error probability. 

I. Introduction 

It is widely accepted that special functions constitute invaluable mathematical tools in almost 
all areas of natural sciences and engineering. In the wide field of digital communications they 
have been extensively used in numerous analytic studies where they typically constitute the 
derivation of explicit expressions for important performance measures possible. Such measures 
correspond to the output signal-to-noise-ratio (SNR), channel capacity, error probability and 
higher-order statistics, to name a few. In addition, the usefulness of special functions has 
became ever more evident over the last decades due to the fact that the majority of them is 
included as built-in functions in popular scientific software packages such as Maple, Matlab 
and Mathematica. As a result, their traditionally laborious computational realisation as well as 
their algebraic complexity have been significantly simplified. 

Among others, the Rice Je-function and the incomplete Lipschitz-Hankel integrals (ILHIs) 
have been shown to be useful in telecommunications since they have appeared in various analytic 
solutions of problems related to wireless communication systems. They were both proposed a 
few decades ago and they are denoted as Ie{k,x) and Zem,nia, z), respectively Q. In more 
details, the Ie{k,x) function is typically defined by a lower incomplete integral which involves 
an exponential function and a modified Bessel function of the first kind and zero order Q . 
Alternative algebraic representations include two series expressions which were reported in 
lf3l . These series are infinite and are represented in terms of the modified Struve function and 
the modified Bessel function of the first kind, lH], 01. Moreover, the author in [[5]| reported 
useful closed-form identities between the Je-function and the Marcum Q-function |[6l-|[9l. The 
usefulness of Ie{k,x) function in digital communications is based on the fact that it has been 
employed in the study of zero crossing rates, in the analysis of angle modulation systems and 
in radar pulse detection and error rates in differentially encoded systems, 0,0, ifTOl . ifTTI . 

In the same context, the ILHIs belong to a class of incomplete cylindrical functions which are 
one of the J^{x), Iv{x), Y^{x), K,y{x) and H,y{x) Bessel functions, namely: the Bessel function 
of the first kind, the modified Bessel function of the first kind, the Bessel function of the second 
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kind, the modified Bessel funetion of the second kind and the Hankel functions, respectively 
[fn , [|4l . In the wide field of electrical engineering, the ILHIs have been largely encountered 
in analytic solutions of numerous problems related to the theory of electromagnetics, while in 
communication theory they have been used in recent investigations associated with the error 
rate analysis of Multiple-Input-Mutiple-Output (MIMO) systems under imperfect channel state 
information (CSI) employing adaptive modulation, transmit beamforming and maximal ratio 
combining (MRC), [IISI-O (and the references therein). 

Nevertheless, even thought the above functions have been proved particularly useful in digital 
communications, it is noticed that they are both neither available in tabulated form, nor are they 
included as built-in functions in the aforementioned popular scientific software packages. As 
a consequence, their corresponding analytical tractability as well as computational realisation 
appear to be rather laborious and cumbersome. Motivated by this, the aim of this work is the 
derivation of analytic expressions for the Rice Ie{k,x) function and the Iem,n{o,,z) integrals. 
Specifically, a remarkably accurate closed-form approximation is firstly derived for the Ie{k,x) 
function. Subsequently, an exact closed-form expression is derived for the Iem,n{a,z) integrals 
for the special case that n is an odd multiple of 1/2, i.e. n + 0.5 G N and a polynomial 
approximation which, as in the case of Ie{k,x) function, is valid for all parametric values. 
Furthermore, upper bounds are derived for the corresponding truncation error of the proposed 
polynomial approximations. Importantly, these bounds are expressed in closed-form and are 
shown to be particularly tight. It is also noted that the derived expressions have a tractable 
algebraic representation which ultimately renders them easily computable and useful in various 
analytical studies in wireless communications. Indicatively, such studies include the derivation 
of explicit expressions for vital performance measures, such as channel capacity and probability 
of error, in the wide field of digital communications over fading channels and the information- 
theoretic analysis of MIMO systems. 

The remainder of this paper is structured as follows: Section II revisits the definition and 
basic identities of the Rice le function and the integrals. Subsequently, Sections III 

and IV are devoted to the derivation of novel analytic expressions for each of these functions, 
respectively. Finally, discussion on the potential applicability of the offered relationships in 
wireless communications along with closing remarks, are provided in Section V. 


4 


II. Definitions and existing representations. 

A. The Rice I e-function 

Definition 1. For k,x E M"*" and 0 < k < 1, the Rice I e-function is defined as, 

pX 

Ie{k,x)= / e~^Io{kt)dt, 

Jo 

where /o(.) is the modified Bessel function of the first kind and zero order, |[T1, [|3ll-[|5l. 


( 1 ) 


An equivalent integral representation to © was given in [[3]|, namely, 


Ie{k, x) = 


1 1 ^—x{l—cos6) 


yjl - k? TT Jo 1 - kcosO 
along with the following alternative series expressions, 


dO, 


( 2 ) 




2Vl - P ^ n\ 2^\/l - k"^ 

^ n=0 




L^+i (xVl - + F_1 (xVl - k^^ 


(3) 


and 


Ie{k,x) = xe 

n=0 

r— OO 


X (1 — k^) 

2fc 

X (1 — 


2 ^ \ n+1 j 
J- 1 


n+1 


{kx) 


r (n +1) 


+ 


n=0 


^n+1 


(fcx) 


2k ; r(n + |) 


(4) 


The notations T„(.) and r(.) denote the modified Struve function and the gamma function. 


respectively ©, [0. It is recalled that Q converges relatively quickly for the case that xylFF 
is large and kx is small, whereas (jH) converges quickly when x\/l — k‘^ is small and kx 
is large. Nevertheless, computing the Ie{k,x) function using these expressions is inefficient 
due to the following three reasons: i) two relationships are required; ii) the above series are 
relatively unstable due to their infinite form; iii) the L„(.) function is not built-in in widely used 
mathematical software packages. 

An adequate way of resolving this issue was reported in [I5]| where the Ie{k,x) function is 
related to the known Marcum Q-function of the first order, namely. 
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Ie{k, x) = ^ [2Q(a, b) - e ""loikx) - l] 

and 

Ie{k, x) = ^ [Q{a, b) - Q{b, a)] 

vi — 

where, 

a = -s/x^J 1 + Vl — 


b = y/xy 1 — \/l — /c^ 


and 


Qi{a,b) = Q{a,b) = xe ^ Io{ax)dx 


whieh was firstly proposed in [0. 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


B. The Incomplete Lipschitz-Hankel Integrals. 

Definition 2. For z G and ^ a,m,n G M, the generalized Incomplete Lipschitz-Hankel 
Integral (ILHI) is defined as, 


Zem,nia,z)= / '"^Zn{x)dx 


( 10 ) 


where Zn{x) denotes one of the cylindrical function^ Jn{x), In{x), Yfix), Kn{x), Hl{x) or 
Hl{x), JU, 0. 


Aeeording to ifT^ . llT3l the parameters m, n, a, z can be also eomplex. An alternative rep¬ 
resentation for the ILHIs of the first-kind modified Bessel funetions in terms of the Marcum 
Q-function was recently reported in IfT^ . namely, 

*The present analysis is limited to the consideration of the In{x) function. 
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m n+1 

i=0 j=0 


+ ^m,n(«)Ql 




+ \/c 


^/x\J a + Va^ - ij 


( 11 ) 


where the eoeffieients y4|„ „(a) and ean be obtained reeursively. As already mentioned, 

the above expression was derived in the eontext of the error rate analysis of MIMO systems 
with imperfect channel state information (CSI). 

III. Novel Expressions For the Rice /c-Function. 

This Section is devoted to the derivation of a tight polynomial approximation for the Rice 
Je-function. 

Proposition 1. For x, fc G M’*' and 0 < k < 1, the following expression holds, 


Ie{k, x) ~ 


r(L + /)Li-2'A;2^7(l + 2/,x) 


( 12 ) 


/!r(L-/ + i)r(/ + i)22' 

where r(x) and 7(0, x) denote the complete and incomplete gamma functions, respectively. 

Proof. It is firstly recalled that a simple and remarkably tight approximation for the In{x 
function was reported in [fT5l eq. (20)], namely. 


In{x) ^ 


T{L + l)L 


1-21 


X\n+2l 


(13) 


^/!r(L-/ + l)r(n + / + l) V2y 

This expression differs from the corresponding Taylor series representation by the following three 
terms: {L +1)/T{L — I + 1) and has been shown to be more stable and better convergent. 

Therefore, by setting n = 0 and x = kt in (11) and substituting in ([1]), it follows that: 


Ie{k, x) ~ 


r(L +/)Li-2'A;2i 


1=0 


/!r(L-/ + i)r(/ + i)22' 


^dt. 


(14) 


Notably, the above integral can be expressed in closed-form in terms of the lower incomplete 
gamma function according to (H eq. (3.381.3)]. As a result, equation (10) is deduced and hence, 
the proof is completed. ■ 
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Remark. As L —>■ oo, equation (10) degenerates into the following exact infinite series rep¬ 
resentation, 


Ie{k, x) 


OO 


E 

z=o 


A:^^ 7(1 + 21, x) 

/!r(/ + 1)22' 


(15) 


which to the best of the authors’ knowledge has not been previously reported in the open technical 
literature. 



Fig. 1. Behaviour of the polynomial approximation in (10) w.r.t k for different values of x and truncation after 20 terms. 


A. A Closed-form Upper Bound for the Truncation Error. 

As already mentioned, equation (10) is an accurate polynomial approximation to the Ie{k,x) 
function. Evidently, the corresponding accuracy increases with L and as with all cases involving 
non-finite series expansions, a closed-form bound for the truncation error is undoubtedly neces¬ 
sary. 


Lemma 1. For G and 0 < k < 1, the following inequality holds, 
















X 


Fig. 2. Behaviour of the polynomial approximation in (10) w.r.t x for different k values and truncation after 20 terms. 


e* < 1 + 



er/(v^\/l - k) 
\/l — k 


erf{y/x\/l + k) 
\/l + k 


L 

/=o 


r(L + /)Li-2'P'7(l + 2/,a;) 
/!r(L-/ + l)r(/ + l)22« 
(16) 


where erf(x) denotes the error function. 


Proof. For the case that (10) is truncated at the term, the corresponding truncation error is 
expressed as 


Q — 


OO 


l=L+l 


r(L + /)Li-2'P'7(l + 2/,x) 
/!r(L-/ + i)r(/ + i)22« 


(17) 


Since as already shown for L —)■ oo, equation (fT^ yields an exact infinite series for the le- 


function, it follows straightforwardly that 


= 


OO 


E 

1=0 


fc^'7(l + 21, x) 

/!!(/ +1)22' 


which can be equivalently written as 


L 


E 

1=0 


r(L + /)Li-2'A;2'7(l + 2/,a;) 
/!r(L-/ + l)r(/ + l)22' 


(18) 























(19) 


et = Ie{k,x) 


L 


E 

1=0 


r(L + /)Li-2^F7(l + 2/,a;) 
/!r(L-/ + l)r(/ + l)22« 


Importantly, the above expression ean be upper bounded by using the tight elosed-form upper 
bound for the Ie{k,x) funetion reported in |[I3 eq. (13)], namely. 


Ie{k,x) < 1 — e ^Io{kx) + 


erf[c^) _ erf{dy/x) 
c d 


( 20 ) 


where c = \/l — /c and d = \/l + k. Henee, by substituting (18) into (17) yields (14), whieh 
eompletes the proof. ■ 


Remark. By setting L ^ cxd in the the finite series in (14) yields a similar elosed-form upper 
bound ean be also obtained for the exaet infinite series in (13), namely. 


e* < 1 + 



er/(v^\/l - k) 
\/l — k 


er/(v^-\/l -f k) 
\/l + k 


Io{kx)e ^ 


^ k^h{l + 2l,x) 

^ /!!(/ + 1 ) 22 ' • 1 1 


B. Numerical Results 

This Seetion is devoted to the validation and analysis of the performanee of the offered results. 
To this end. Figure 1 illustrates the behaviour of Ie{k,x) funetion as funetion of k for different 
values of x and with the series truneated after 20 terms. One ean see the monotonieally inereasing 
behaviour of the funetion with respeet to k. Notably, the series eonvergenees partieularly quiekly 
sinee for all seenarios with the eorresponding truneation after 20 terms, the both the involved 
absolute error =| Ie{k, x)—Ie{k, x) \ and the absolute relative error =| Ie{k, x)—Ie{k, x) \ 
/Ie{k,x) is signifieantly less than 10“^°. 

Figure 2 depiets the behaviour of (10) as a funetion of x for different values of k and for 
truneation after 20 terms. In the small value regime, the differenee between the three seenarios 
is rather small; however, as x inereases the monotonieally inereasing behaviour of the funetion 
beeomes more evident. Furthermore, the involved absolute and absolute relative errors are also 
very low, e^, ear < 10“®, in every possible seenario. 
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Fig. 3. Behaviour of the closed-form expression in (20) w.r.t m for z = A, a = 1.8 and different values of n. 


IV. Novel Expressions For the ILHIs. 

This Section is devoted to the derivation of a closed-form expression for the case that n is 
a half integer as well as a tight polynomial approximation which is valid for all values of the 
involved parameters. 

A. A Closed-Form Expression for the case that n + | G R 


Theorem 1. For a, m, z G M, a > 1 and n + ^ G N, the following expression holds, 


I Cr, 


[a, z = 




{n + k-l)\ 


^ y/iikl (n - k - 1) !2^+ 


where A = m — k + 


iMc7(A(a-1)^) , ^ ^^^+17(71, (a + l)^)l 

(0-1).^ “imp—/ 

( 22 ) 


Proof. It is recalled that for the case that n + 0.5 G M, the modified Bessel function of the 
first kind is expressed in closed-form by dH eq. (8.467)], namely, 
















II 




liX = 




{n + k)\ [(-1)V + (-l)’"+ie-^] 


(23) 


k=o v^fc!(n -/c)!(2a;)^+2 

To this effect, by making the necessary change of variables and substituting in (8), one obtains, 




2 


^ 2 


[n + k - l)\2 

^ 0Ffc! {n-k-l)\ 


- 1 ) 


k I ^A^-ax„x 


e^dx + {—lY^'^ I x''e 

Jo 




(24) 


Notably, the integrals in (22) belong to the family of gamma special functions. Therefore, after 
some basic algebraic manipulation and with the aid of flU eq. (3.381.3)], equation (20) is deduced 
and therefore, the proof is completed. ■ 


B. A Novel Series Representation 


Proposition 2. For m,a,n,z G M, the following expression holds, 


, \ r(Z7 + Z)T^ ^*7(m + 77. + 2/+ 1, az) 

f e™,„(a, ^ ^ ^ ^ ^ l^2^+2l^ra+n+2l+l ^ 

where r(a;) and 7(0, x) denote the complete and incomplete gamma functions, respectively. 


Proof. Recalling the Gross’ polynomial approximation for the /^(x) function, substituting (11) 
in (8), recalling that a! = r(a + 1) and performing some basic algebraic manipulations, one 
obtains 


k ^ ^ 


r(L + /)L 


l—2ln—n—2l 


1=0 


/!(L-/)!r(n + / + l)a™+^+2^+i Jo 


^m+n+2l ^-ax 


(26) 


By making the necessary change of variables and using [HI eq. (3. 381.3)] yields (23) and thus, 
the proof is completed. ■ 


Remark. By recalling that (11) becomes exact as L —)■ 00 , it follows straightforwardly that 
can be expressed by the following explicit expression. 
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Fig. 4. Behaviour of the closed-form expression in (21) w.r.t z for a = 1.8, m = 1.2 and different values of n. 


1=0 


7 (m + n + 2/ + 1, az) 
nT{n + I + i)2^+2ia”^+^+2«+i' 


To the best of the authors’ knowledge the above expressions are novel. 


(27) 


C. A Closed-Form Bound for the Truncation Error 

A closed-form bound for the corresponding truncation error can be derived by following the 
same methodology as in Lemma 1. 


Lemma 2. For m, a,n, z G R, the following inequality holds, 


et < 


Ln-ij + i 


k=0 


{k + [n - + l)\ 

y/nk\ — + \ — k)\2^^^ 


(a — 1)-^ (a + 1)-^ 


L 


1=0 


r(L + l)L^ ^^ 7 (m + n + 2 Z + 1, az) 
1\{L - l)\T{n + I + i)2»^+2;a"^+^+2'+i 


(28) 
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Fig. 5. Behaviour of the closed-form expression in (23) w.r.t m for a = 2.2, 2 = 5 and different values of n for truncation 
after 700 terms. 


where A = m — k + 


Proof. By truncating the series in (23) after L terms, the eorresponding truneation error is 
given by 

\ r(ii + k)L} + fi + 2/ + 1, az) 

^ 1 \{L - /)!r(n + I + i)2^+2ia^+n+2i+i ^ ’ 

Importantly, given that 

oo oo L 

/(^) = 

i=L+l 1=0 I 

and based on the derived series in (25), equation (27) is equivalently re-written as 


e* = 


L 

k • 2 ) ^ ^ 

1=0 


T{L + l)L^ ^'7(777 -I- n -f 2 / -f 1 , az) 
1 \{L - /)!r(n + I + l)2^+2lam+n+2l+l 


(31) 
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It is recalled here that the Iem,n{cL,z) integral is monotonically decreasing with respect to n. 
Thus, by assuming an arbitrary positive real scalar b G the following inequality holds, 


and thus Je^ „_,_i (a, < Iem,n{a,z). Therefore, by recalling the closed-form representation in 

Theorem 1 and making use of the floor function principle in order to formulate the identity 
n> \n — 0.5J + 0.5, the /em,n(a, z) function can be upper bounded Vn e M+ as follows. 


{k+[n-l\ + l)\ 


lerania, z) < > - 

’ to V^kl {[n - l\ + I - 

/ (g- 1)^) . 1^+11 7 (A (a + 1)^) 

^ ^ ’ (a + 1)^ 


X 


(33) 


where A = m—k+^. Therefore, by inserting (30) into (28) yields (26) and therefore completes 
the proof. ■ 


Remark. By recalling that (23) becomes exact as L —)■ cxo, a corresponding closed-form bound 
for the truncation error of (25) is straightforwardly deduced, namely. 


e* < 


VnA\ 


{k+ [n- |J)! 


^ 0rfc! ([n - 2j - fc)!2^+2 


^ ^ ^ (g + l)^) 

(a — 1)-^ (a -f 1)-^ 


L 


-E 

1=0 


7 (m + n + 21 + az) 
/!r(n + I + i)2^+2«a™+^+2'+i 


(34) 


where again A = m — k + ^. 

It is recalled that the offered results can be used in numerous areas in natural sciences and 
engineering including wireless communications with emphasis in digital communications over 
fading channels llUTl - Gbll and the references therein. 
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Fig. 6. Behaviour of the polynomial approximation in (23) w.r.t 2 for a = 1.4, m = 2.2 and different values of n for truncation 
after 700 terms. 


D. Numerical Results 

The validity of the derived elosed-form expression and the general behaviour of the offered 
bounds are illustrated in Figures 3 — 6. Figure 3 depiets the behaviour of the elosed-form 
expression in (20) with respect to m for different half integer values of n and assuming that 
z = 4: and a = 1.8. One can observe and verify the monotonically decreasing behaviour of 
the function with respect to n. This is also observed in Figure 4 where the behaviour of (20) 
is depicted with respect to 2 ; for a = 1.8 and m = 1.2. In the same context, Figures 5 and 6 
illustrate the behaviour of the in (23) for the case of arbitrary values of n and for 

truncation after 700 terms. The former is represented w.r.t m for a = 2.2 and z = 5 while the 
latter is plotted w.r.t 2 ; for a = 1.2 and m = 2.2. The usefulness of the offered expressions 
is evident since (20) is an exact closed-form expression while (23) is simple to use and very 
accurate. It is additionally useful for its capability to account for the sensitivity of n. 
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V. Conclusion 

In this work, novel analytic results were presented for the Rice Je-function and the Incomplete 
Lipschitz Hankel Integrals of the modified Bessel function of the first kind. The offered results 
have a convenient algebraic representation and therefore they are easy to handle both analytically 
and numerically. This is sufficiently advantageous since it constitutes them suitable for utilization 
in various studies related to the analytical performance evaluation of of digital communications 
over fading channels, among others. 
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